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Using a sample of (225.2 ± 2.8) x 10 6 J/ip events collected with the BESIII detector, we present 
results of a study of J/ip — > jAA and report the first observation of the isospin violating decay 
J/ip —¥ AS + c.c, in which X" decays to jA. The measured branching fractions are B{J/ip — > AS ) 
= (1.46 ± 0.11 ± 0.12) x 10~ 5 and B(J/ip AS ) = (1.37 ± 0.12 ± 0.11) x 10" 5 . We search for 
/1(1520) —t jA decay, and find no evident signal, and an upper limit for the product branching 
fraction B(J/ip ->■ /M(1520) + c.c.) x S(yl(1520) -> jA) < 4.1 x 10~ 6 is set at the 90% confidence 
level. We also report the observation of n c — > AA in J/tp — >• yn c , n c — > AA and measure the 
branching fraction B(rj c -> AA) = (1.16 ± 0.12(stat) ± 0.19(syst) ± 0.28(PDG)) x 10~ 3 . 

PACS numbers: 13.25.Gv, 12.38. Qk, 13.60.Rj, 14.20.Jn 



I. INTRODUCTION 

The study of charmonium meson decays into baryon 
pairs is an important field that intersects particle and 
nuclear physics, and provides a novel means for ex- 
ploring various properties of baryons The decay 
J/ip — > AS + c.c. is an isospin symmetry breaking de- 
cay, and a measurement of its branching fraction will help 
elucidate isospin-breaking mechanisms in J/ip — ¥ B^Bg 
decays 0, S|. Until now, only an upper limit on the 
branching fraction of B(J/ip -> AS + c.c.) < 1.5 x 1CT 4 
has been set at the 90% confidence level (C.L.) by the 
MarkI Collaboration, based on a study of J/tp — > ^AA 

The electromagnetic decays of hyperons A* — > jA pro- 
vide clean probes for examining the internal structure of 
A* hyperon resonances [5]. For example, predictions for 
the radiative decay 71(1520) — > jA have been made in a 
number of frameworks including: a nonrelativistic quark 
model [1, 0] ; a relativistic constituent quark model [H ; 



the MIT bag model ; the chiral bag model Q ; an alge- 
braic model of hadron structure [loj ; and a chiral quark 
model [Tlj |. I n contrast, experimental measurements have 
been sparse (l2lJl5l |. The radiative decays A* — > jA can 
be studied with J/ip — > jAA events. 



The J/ip — > jAA events can also originate from radia- 
tive J/ip — > 77y c decays followed by rj c decays to AA. To 
date, rj c — > AA has only been observed in B ± — > AAK^ 
decays by the Belle experiment flo| . A measurement of 
r/ c — > A A in J/ip radiative decays provides useful infor- 
mation in addition to Belle's measurement in B decays. 



In this paper, we report the first observation of the 
isospin violating decay J/ip — > AS + c.c, a new mea- 
surement of the branching fraction for rj c — >• AA and the 
results of a search for the radiative decay A(1520) — > jA. 
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II. DETECTOR AND MONTE CARLO 
SIMULATIONS 

The analysis is based on analyses of J/ip — > jAA 
events contained in a sample of (225.2 ± 2.8) x 10 6 J/ip 
events [l?} accumulated with the Beijing Spectrometer 
III (BESIII) operating at the Beijing Electron-Position 
Collider II (BEPCII) [lj]. 

BEPCII is a double ring e + e~ collider with a design 
peak luminosity of 10 33 cm~ 2 s _1 with beam currents of 
0.93 A. The BESIII detector consists of a cylindrical core 
comprised of a helium-based main drift chamber (MDC), 
a plastic scintillator time-of-fiight (TOF) system, and 
a CsI(Tl) electromagnetic calorimeter (EMC) that are 
all enclosed in a superconducting solenoidal magnet that 
provides a 1.0 T axial magnetic field. The solenoid is 
supported by an octagonal flux-return yoke that con- 
tains resistive-plate-chamber muon-identifier modules in- 
terleaved with plates of steel. The acceptance for charged 
particles and photons is 93% of 47r sr, and the charged- 
particle momentum and photon energy resolutions at 1 
GeV are 0.5% and 2.5%, respectively. 

The responses of the BESIII detector are modeled with 
a Monte Carlo (MC) simulation based on GEANT4 0, 
l20l |. evtgen [2]| is used to generate J/ip -> AE° + c.c. 
events with an angular distribution of 1 -I- a cos 2 9, where 
is the polar angle of the baryon in the J/ip rest frame 
and a is a parameter extracted in fits to data described 
below. The J/ip — > jr/c decays are generated with an 
angular distribution of 1 + cos 2 9 1 and a phase-space dis- 
tribution for rj c — > A A, and effect of spin-correlation is 
not considered in the MC simulation for r\ c — > AA de- 
cay. Inclusive J / ip decays are produced by the MC event 
generator kkmc [23], the known J/ip decay modes are 
generated by evtgen [2l[ with branching fractions set 
at their Particle Data Group (PDG) world average val- 
ues [23| . and the remaining unknown decays are gener- 
ated with LUNDCHARM [2J|. 



III. DATA ANALYSIS 

Charged tracks in the BESIII detector are recon- 
structed from track-induced signals in the MDC. We se- 
lect tracks within ±20 cm of the interaction point in the 
beam direction and within 10 cm in the plane perpendic- 
ular to the beam; the track directions are required to be 
within the MDC fiducial volume, |cos#| < 0.93. Candi- 
date events are required to have four charged tracks with 
net charge zero. The A A pair is reconstructed using the 
A — > pn~ , and A — > pir + decay modes. We loop over all 
the combinations of positive and negative charged track 
pairs and require that at least one (ptt~)(Ptt + ) track hy- 
pothesis successfully passes the A and A's vertex finding 
algorithm. 

If there is more than one accepted (p7r~)(p7r + ) combi- 
nation in an event, the candidate with minimum value of 
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FIG. 1: A scatter plot of M pn - versus M p7T + for selected 
candidate events. 



(Mpn+ ) and Ma {Ma) are the measured mass and its ex- 
pected value. Since there are differences in the detection 
efficiencies between data and the MC simulation for low- 
momentum proton and antiprotons pEj , we reject events 
containing any proton or antiproton track candidate with 
momentum below 0.3 GeV/c. 

Electromagnetic showers are reconstructed from clus- 
ters of energy deposits in the EMC. The energy deposited 
in nearby TOF counters is added to improve the re- 
construction efficiency and energy resolution. Showers 
identified as photon candidates are required to satisfy 
fiducial and shower- quality requirements: e.g., showers 
in the barrel region (| cos f?| < 0.80) must have a min- 
imum energy of 25 MeV, while those from end caps 
(0.86 < |cos0| < 0.92) must have at least 50 MeV. To 
suppress showers generated by charged particles, we re- 
quire that the photon candidate direction is at least 5° 
away from its nearest proton and charged pion tracks, 
and at least 30° away from the nearest antiproton track, 
since more EMC showers tend to be found near the di- 
rection of the antiproton. This requirement decreases 
the signal efficiency by 18% for J/ip -t AS°(S° -> -fA) 
compared to that for J/ip — > AS Q (U — > jA) since the 
photon from the radiative S° — > 7/I decay is closer to 
the direction of the antiproton. Requirements on the 
EMC cluster timing are used to suppress electronic noise 
and energy deposits that are unrelated to the event. 
A four-constraint (4C) energy-momentum conservation 
kinematic fit is performed to the ^AA hypothesis. For 
events with more than one photon candidate, the combi- 
nation with the minimum \\ c ^ s selected. In addition, we 
also require xlc < 45 in order to suppress backgrounds 
from the decays J/ip -» AA, and AAtt . 

A scatter plot of M p7r ~ versus M pw + for events that 
survive the above requirements is shown in Fig.[T|), where 
a cluster of A and A signals is evident. To select J/ip — > 
j A A signal events, we require \M pn - — M A \ < 5 MeV/c 2 
and \M p% + - M A \ < 5 MeV/c 2 . An M 2 (jA) (vertical) 
versus M 2 (7^1) (horizontal) Dalitz plot for these events 
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FIG. 2: (a) An M 2 ('yA) (vertical) versus M 2 (-yA) (hori- 
zontal) Dalitz plot for selected events and the (b) 7/I , (c) 
jA and (d) AA invariant mass distributions for the selected 
J/ip — > jAA event sample. 




1.2 1.25 
My A (GeV/c 2 ) 




M (GeV/c") 



is shown in Fig. [2] (a); the 7^1 and mass spectra are 
shown in Fig. [5] (b) and (c). Prominent signals of the 
S° and S a , corresponding to J/ip — ► AS + c.c. decays, 
are observed. On the other hand, no obvious signal for 
/1(1520) — > j A is seen. A clear r\ c signal can be seen 
in the AA mass spectrum shown in Fig. [5] (d), while no 
significant enhancement at other A A masses is evident. 

For the J/ip — > AS + c.c. study, we apply the same 
requirements to a sample of 225 x 10 6 MC-simulated 
inclusive J/ip events and find that the primary back- 
grounds come from J/ip — > AA, S°S° and AAir de- 
cays, where either a cluster in the EMC unrelated to the 
event is misidentified as a photon candidate or one of 
the photons from the E°S° or ir° decay is undetected in 
the EMC. Normalized M(-yA) and M(-yA) distributions 
from the events that survive the application of the 4C 
kinematic fit, shown as dotted and dashed histograms 
in Figs. [3] (a) and (b), show no sign of peaking in the 
S° or S° mass regions. Another potential source of 
background is from J /ip — ¥ jr/c (rj c — > AA) decay and 
nonresonant J/ip — > "fAA, which contribute a smooth 
background under the signal region, shown as dot-dashed 
curves in Figs. [3] (a) and (b). The expected backgrounds 
are 105 ± 10 (95_± 9) events in the S° (S°) signal re- 
gion for J/ip ~> AS (J '/ip AS ) as listed in Table H 
The signal region is defined as being within ±3<r of the 
nominal £° (S°) mass. It should be noted that the back- 
ground events from the nonresonant J/yJ —¥ ^AA are not 
counted and are accounted for by the floating polynomial 
function discussed below. 

Unbinned maximum likelihood (ML) fits are used to 
determine the AS and AS event yields. The signal 



FIG. 3: The results of the fit for the E° (a) and E° (b). 
The points with error bars are data. The fit results are 
shown by the black solid curves. The light (red) solid curves 
are the signal shapes. The (blue) dotted histograms are 
from the normalized J/ip — > AA background; the (green) 
dashed histograms are from the normalized Jip — * 
background. The (magenta) dot-dashed curves show the non- 
resonant (phase-space) background polynomial. 



probability density function (PDF) for S° (S°) from 
J/ip — > AS (AS ) is represented by a double-Gaussian 
function with parameters determined from the MC simu- 
lation except for the Gaussian widths, which are allowed 
to float. Backgrounds from J/ip AA and £°XJ a are 
fixed to their MC simulations at their expected intensi- 
ties. The remaining background is described by a second- 
order polynomial function with parameters that are al- 
lowed to float. The fitting ranges for both the £° and the 
£° are 1.165- 1.30 GeV/c 2 . Figures [3] (a) and (b) show 
the results of the fits to S° and S°. The fitted yields 
are 308_± 24 and 234 ± 21 signal events for J/tjj AS 
and AS , respectively. The goodness of fit is estimated 
by using a % 2 test method with the data distributions 
regrouped to ensure that each bin contains more than 
10 events. The test gives x 2 /n.d.f = 28.1/37 = 0.76 
for J/i/j -> AS and x 2 /n.d./ = 43.5/37 = 1.2 for 
J/ip — > AS , where n.d.f. is the number of degrees of 
freedom. 

In the higher jA (jA) invariant mass regions, shown in 
Figs. H](b) and (c), no obvious signals for yl(1520) — > jA 
(/1(1520) — > jA) are evident. We require that the in- 
variant mass of AA is less than 2.9 GeV/c 2 to further 
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FIG. 4: The results of the fit for the 71(1520). The points 
with error bars are data. The fit result is shown by the black 
solid curve; the (magenta) dashed curve is the background 
polynomial and the (red) light solid curve is the 71(1520) signal 
shape. [Here the M(yA) and M(jA) mass distributions are 
combined.! 



FIG. 5: The r\ c mass distribution and fit results. Points with 
error bars are data. The fit result is shown as a black solid 
curve, the (red) light solid curve is the signal shape, the (blue) 
dashed curve is the combined incoherent background from the 
J/ip — > E°E° , AE° + c.c, the (magenta) dot-dashed-curve is 
the nonresonant background. 



suppress combinatorial backgrounds from J/tp — > AA, 
E°E , AAtt and J/ip ->■ 7r/ c (r/ c AA) decays. Af- 
ter the above requirement, only 14 ± 1 events from 
these background decay modes remain. In the surviving 
combined M{jA) and M(-yA) mass spectrum, shown in 
Fig. El there is no evidence for a /1(1520) signal above ex- 
pectations for a phase-space distribution of J/ip — > yAA. 

In the ML fit to the Fig. H distribution, the /1(1520) 
signal PDF is represented by a Breit-Wigner (BW) func- 
tion convolved with a double-Gaussian resolution func- 
tion, with parameters determined from the fit to the S° 
data. The shape for the nonresonant background is de- 
scribed by a second-order polynomial function, and the 
background yield and its PDF parameters are allowed to 
float in the fit. The mass range used for the A(1520) fit is 
1.35-1.70 GeV/c 2 . Figure H shows the result of the fit to 
yl(1520), which returns a /L(1520) signal yield of 31 ± 24 
events. The goodness of fit is x 2 /n.d.f = 45.9/45 = 1.02. 
Using a Bayesian method, an upper limit for the number 
of /1(1520) signal events is determined to be 62.5 at the 
90% confidence level (C.L.). The signal yields and the 
efficiencies for the analyses of J/ip — > A2J° (AS ) and 
ylyl(1520) + c.c. are summarized in Table |U 

For the J/ip — > ^r\ c (j]c AA) analysis, the dom- 
inant backgrounds remaining after event selection are 
from J/ip — > U°IJ and AZ!° + c.c The expected num- 
ber of events in the signal region from these two sources 
is 637 ± 52, as listed in Table Q] These backgrounds are 
incoherent (i.e., do not interfere with the rj c signal am- 
plitude). In addition, there is an irreducible background 
from nonresonant J/ip — > yAA decays that is potentially 
coherent with the signal process (i.e., may interfere with 
the rj c signal amplitude). 

For the rj c fit, the combined incoherent background 
is fixed to the shape and level of the MC simulation. 
The PDF for the coherent nonresonant background is 



TABLE I: For each decay mode, the number of signal events 
(Ns), the number of expected background events (Nb) in the 
signal region (nonresonant J/tp —¥ yAA background is ex- 
cluded), and the MC efficiency (e) for signal are given. The 
error on Ns is statistical only, and the signal regions are de- 
fined to be within ±3cr of the nominal X" and 71(1520) masses. 



Modes 


N S 


N B 


e(%) 




AS°(E° -> 7 /l) 




308 ± 24 


105 ± 10 


21.7 




AE°(E° -» j A) 




234 ± 21 


95 ± 9 


17.6 




/L4(1520) + c.c.(A(1520) - 


•> yA) 


31 ± 24 


14 ± 1 


18.8 




■yr] c (r] c -5- AA) 




360 ± 38 


637 ± 52 


19.8 



described by a second-order polynomial, with yield and 
shape parameters that are floated in the fit. For the line- 
shape for rj c mesons produced via the Ml transition, we 
use (E^ x BW(rn) x damping(£J 7 )) <g> Gauss(0, a), where 

— M 2 —™ 2 

m is the AA invariant mass, E*y 



2M 



J/<l> 



is the en- 



ergy of the transition photon in the rest frame of J/ip, 
damping(.E 7 ) is a function that damps the divergent low- 
mass tail produced by the factor, and Gauss(0, a) is a 
Gaussian function that describes the detector resolution. 
The damping function used by the KEDR [26{ collabo- 
ration for a related process has the form 



EqE^ + (Eq 



(1) 



Mi 



where Eq — — 'J' lc is the peak energy of the transi- 
tion photon. On the other hand, the CLEO experiment 
damped the Ejt term by a factor exp(-E^/8(3 2 ), with 
P — 65 MeV [27| , to account for the difference in overlap 
of the ground state wave functions. We use the KEDR 
function in our default fit and use the CLEO function 
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as an alternative. The difference between the results ob- 
tained with the two damping functions is considered as a 
systematic error associated with uncertainties in the line 
shape. In the fit, the mass and width of r\ c are fixed to 
the recent BESIII measurements: M(r]^ = 2984.3 ± 0.8 
MeV/c 2 and r(ry c ) = 32.0 ± 1.6 MeV [11], and interfer- 
ence between the nonresonant background and the rj c res- 
onance amplitude is neglected [2!| . The mass range used 
for the rj c fit is 2.76 — 3.06 GeV/c 2 . Figure [5] shows the re- 
sult of the fit to rjc, which yields (360 ± 38) signal events. 
The goodness of the fit is x 2 /n.d.f = 42.7/43 = 0.99. 
The signal yield and efficiency are summarized in Table HI 



IV. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties on the branching fraction 
measurements are summarized in Table [Til The system- 
atic uncertainty due to the charged tracking efficiency has 
been studied with control samples of J/ip —¥ pK~ A + c.c. 
and J/ip — > AA decays. The difference in the charged 
tracking efficiency between data and the MC simulation 
is 1% per track. The uncertainty due to the A and A 
vertex fit is determined to be 1% for each A by using 
the same control samples. The uncertainty due to the 
photon reconstruction is determined to be 1% for each 
photon 



17j . The uncertainties due to the kinematic fit 



are determined by comparing the efficiency as a function 
of xlc v& l ue f° r the MC samples and the control samples 
of J/ip -t AA and J/ip 



S°S° events, in which zero 



and two photons are involved in the final states. The 
differences in the efficiencies between data and MC sim- 
ulation are 2.1% and 2.3% from the studies of J/ip — > A A 
and J/ip — > S°S° events, respectively; we use 2.3% as 
the systematic error due to the kinematic fit. 

The signal shape for the S° (S°) is described by a 
double-Gaussian function and the widths are floated in 
the nominal fit. An alternative fit is performed by fixing 
the signal shape to the MC simulation, and the system- 
atic uncertainty is set based on the change observed in 
the yield. In the fit to .4(1520), since the shape of the 
signal is obtained from MC simulation, the uncertainty 
is estimated by changing the mass and width of 4(1520) 
by 1 standard deviation from their PDG world average 
values [23]. This systematic error is determined in this 
way to be 4.8%. 

In the r] c fit, the mass resolution is fixed to the MC sim- 
ulation; the level of possible discrepancy is determined 
with a smearing Gaussian, for which a nonzero a would 
represent a MC-data difference in the mass resolution. 
The uncertainty associated with a difference determined 
in this way is 1.1%. Changes in the mass and width of 
the rj c used in the fit by 1 standard deviation from the 
recently measured BESIII values [28], produce a relative 
change in the signal yield of 6.4%. As mentioned above, 
damping functions from the KEDR and CLEO collabora- 
tions were used in the fit to suppress the lower mass tail 
produced by the factor; the relative difference in the 



yields between the two fits is 3.9%. The 7.6% quadrature 
sum of these uncertainties is used as the systematical er- 
ror associated with uncertainties in rj c signal line-shape. 

For the measurement of the J/ip -» AS (AS ), the 
expected number of background events from the decays 
of J/ip -> AA and S°S° is fixed in the fit. To estimate 
the associated uncertainty, we vary the number of ex- 
pected background events by 1 standard deviation from 
the PDG branching fraction values (23[, which gives an 
uncertainty of 0.6% (0.4%) for the J/ip -> AS (AS ). 
In the ML fit to r\ c , the incoherent backgrounds from 
J/ip — > S°S° and AS + c.c. are also fixed at their ex- 
pected numbers of events. The uncertainty associated 
with this is determined by changing the number of ex- 
pected incoherent background events by 1 standard de- 
viation of the PDG branching fraction values [23[ for 
the J/ip — > S°S° and the measured value for J/ip — > 
AS + c.c. from the analysis reported here; the resulting 
change in the r\ c signal yield is 12.8%. 

The uncertainty due to the nonresonant background 
shape for each mode has been estimated by changing 
the polynomial order from two to three. The systematic 
uncertainties due to the fitting ranges are evaluated by 
changing them from 1.165- 1.30 GeV/c 2 to 1.165-1.25 
GeV/c 2 (S° and S°), from 1.35- 1.70 GeV/c 2 to 1.38- 
1.67 GeV/c 2 (4(1520)) and from 2.76 - 3.06 GeV/c 2 to 
2.70 - 3.06 GeV/c 2 (rj c ). The changes in yields for these 
variations give systematic uncertainties due to the choices 
of fitting ranges, as shown in Table [TTJ 

The electromagnetic cross sections for AS + c.c. pro- 
duction through direct one-photon exchange and J/ip de- 
cay in e + e~ can be inferred using the factorization hy- 
pothesis to be [30( | 



cr(e + e 



AS ) a(e+e- -> 7* 



o\c c 



J/ip -> AS ) a(e+e~ -> J/ip -> ' 



(2) 



3.097 GeV, a(e 



N 



Neglecting interference between e + e — > 7* 
and e + e~ — > J/ip — > one can obtain, at y/s = 

+e~ -> J/ip -> n + n~) = B(J/ip -> 
168 ±3.2 nb, where Nj/^ and C are the 
number of total J/ip events (225.2 ± 2.8 x 10 6 ) and the 
corresponding integrated luminosity (79631 ± 70(stat.)± 
796(syst.)) nb 17\, respectively. At ^fs = 3.097 GeV, 



.7/ ,/■ 



OBorn(e + e" 



7* — > it [i ) is 9.05 nb. From this we 



estimate the relative ratio of the QED background from 
e+e" -> 7* -> AS + c.c. to be (5.4 ± 0.1)% of our mea- 
sured yield of J/ip — > AS + c.c. events. Therefore, we 
adjust our result be a factor of 0.946 when we determine 
the J/ip — ► AS + c.c. branching fraction value; we use 
0.1% as a systematic error due to the uncertainty in this 
correction factor. 

The angular distribution of the baryon in J/ip — > BsBg 
decay is expected to have a 1 + a cos 2 6 behavior. Fig- 
ures [6] (a) and (b) show the distributions of cos 9 for 
A (J/ip -> AS ) and 4 (J/ip -> AS ), respectively, 
after correcting the signal yields for the detection effi- 
ciency. A simultaneous fit to the angular distributions 
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TABLE II: Summary of systematic errors for the branching fraction measurements (%). 





J/tp -+ AE° 


J/tp -5- AS 


J/ip -> AA(1520) + c.c. -> 7/M 


J/ip — ► 7r? c — >• 7AA 


Photon detection 


1 


1 


1 


1 


Tracking 


4 


4 


4 


4 


A and A vertex fits 


2 


2 


2 


2 


4C kinematic fit 


2.3 


2.3 


2.3 


2.3 


Signal shape 


1.3 


2.6 


4.8 


7.6 


Fitting range 


1.6 


0.9 


1.4 


1.4 


a 


5.5 


5.1 


10.2 




Fixed backgrounds 


0.6 


0.4 




12.8 


Nonresonant background shape 


0.3 


0.1 


1.9 


1.7 


QED correction factor 


0.1 


0.1 






Cited branching fractions 


0.8 


0.8 


0.8 


0.8 


Number of J /if) 


1.3 


1.3 


1.3 


1.3 


Total systematic uncertainty 


8.0 


7.9 


12.6 


16.0 



for A and A returns the value a = 0.38 ± 0.39. The 
detection efficiencies are determined with MC simula- 
tion for J/ip — > AS ± c.c. using a = 0.38 in the sig- 
nal MC generator. To estimate the uncertainty orig- 
inating from the parameter a, we generate MC sam- 
ples for a — 0.38 and for other values in the range 
0.0 ~ 0.77. The maximum difference is 5.1% (5.5%) for 
J/ip — > AS (AS ) and is taken as a systematic error. 
For J/ip — > 7L4(1520)+ c.c. decay, the detection efficiency 
is obtained with a phase-space MC simulation. We gen- 
erate MC samples for a = and a — 1 to estimate the 
uncertainty due to the unknown parameter a. The dif- 
ference of efficiency of 10.2% is taken as systematic error 
for the J/ip -> A4(1520) + c.c. 




FIG. 6: The corrected distributions of cos 9 for A from 
J/ip -»■ AU° decay (a), for A from J/ip -» AS decay (b). 
The curves in (a) and (b) present the fits to the function 
1 + a cos 2 9. The goodness of the fits are \ 2 /n.d.f = 21/18 = 
1.2 for A and x 2 /n.d.f = 29/18 = 1.6 for A. 



The branching fraction for the A — > pir decay is taken 
from the PDG [23|]; the 0.8% uncertainty is taken as a 
systematic uncertainty in our measurements. The uncer- 
tainty in the number of J/ip decays in our data sample 
is 1.3% [l7|. The total systematic uncertainties for the 
branching fraction measurements are obtained by adding 



up the contributions from all the systematic sources in 
quadrature as summarized in Table HT1 

V. RESULTS AND DISCUSSION 

The branching fractions are calculated with B = 
Ns/{Nj/^eBp V ), where N$ and e are the number of sig- 
nal events and the detection efficiency, listed in Table Q] 
Here Nj/^ = (225.2±2.8) x 10 6 [13 is the number of J/ip 
events, and £> p7T is the branching fraction of the A — > pit 
taken from the PDG [23j . The calculated branching frac- 
tions, along with the PDG 23] limits, are listed in Ta- 
ble [TTT. 

TABLE III: Branching fractions (10 ) from this analysis, 
where the first errors are statistical and the second ones are 
systematic, and the PDG values [23| for comparison. The 
upper limits are at the 90% C.L.. 



J/ip decay mode 


BESIII 


PDG 


AE° 




1.46 ±0.11 ±0.12 


< 7.5 


AE° 




1.37 ±0.12 ±0.11 


< 7.5 


~tVc(Vc AA) 




1.98 ±0.21 ±0.32 




/L4(1520) + c.c.(A(1520) - 


^■yA) 


< 0.41 





Our measurement of the branching fraction for J /ip — > 
AS + c.c. decay can shed light on the SU(3) break- 
ing mechanism. The amplitude for J/ip decay to a 
pair of octet baryons can be parametrized in terms of 
a SU(3) singlet A, as well as symmetric and antisym- 
metric charge-breaking (D, F) and mass-breaking (D , 
F ) terms, as described in Refs. 0, H, HU and listed 
in Table IIV[ where 5 is used to designate the rela- 
tive phase between the one-photon and gluon-mediated 
hadronic decay amplitudes. According to these ampli- 
tude parametrizations the J/ip — > AS ± c.c. branching 
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TABLE IV: Amplitude parametrizations from @, HI, [HJ for 
J/tp decay to a pair of octet baryons. General expressions in 
terms of a singlet A, as well as symmetric and antisymmetric 
charge-breaking (D, F) and mass-breaking terms (D , F ) are 
given. Here S is the relative phase between one-photon and 
gluon-mediated hadronic decay amplitudes. Except for the 
branching fraction for J/tp — ¥ AS +c.c. decay (marked with 
an asterisk) from this measurement and for J/tp — > pp, nn 
from the recent BESIII measurements [52] , the other branch- 
ing fractions (B) are taken from the PDG [23| . 



Decay mode 


Amplitude 


B(xlO- a ) 


PP 


A + e iS (D + F) + D' + F' 


(2.112 ± 0.031) [32] 


nn 


A-e iS (2D) + D' +F' 


(2.07 ± 0.17) [321 


£+£- 


A + e iS (D + F) - 2D' 


(1.50 ±0.24) 




A + e lS (D) - 2D' 


(1.29 ± 0.09) 


S°S° 


A - e iS (2D) +D' -F' 


(1.20 ±0.24) 




A + e iS (D - F) + D' - F' 


(0.85 ± 0.16) 


AA 


A - e l5 (D) +2D' 


(1.61 ± 0.15) 


AS (AS ) 


(V3D) 


(0.014 ±0.002)* 



J/tp -> AS , J/tp -> AE° and J/tp -> jr] c (n c ^ AA ) 
are measured for the first time as: B{J/tp — > AU°) = 
(1.46 ± 0.11 ± 0.12) x 10" 5 , B(J/ip -> AS ) = (1.37 ± 
0.12 ± 0.11) x 10" 5 and B{ J/tp -> Th) x B{Vc -> AA) = 
(1.98±0.21±0.32) x 10 -5 , respectively, where the uncer- 
tainties are statistical and systematic. Using the PDG 
value [H[ for J/tp — » r yrj c , we obtain B(r] c — > AA) — 
(1.16 ± 0.12 ± 0.19 ± 0.28 (PDG)) x 10~ 3 , where the 
third error is from the error on B(J/tp — > JT} C )- Us- 
ing B^ — > AAK^ decay the Belle experiment measured 
B(rj c -> AA) = (0.87lg;|ttg;S±0.27 (PDG)) x 10" 3 [H, 
which is consistent with our result within error. No evi- 
dence for the decay of J/tp —> A4(1520) + c.c. is found, 
and an upper limit for the branching fraction is deter- 
mined to be B(J/tp -> AA(1520) + c.c.) x B(vl(1520) -> 
j A) < 4.1 x 10 -6 at the 90% confidence level. Results 
are listed in Table IIIII and compared with previous mea- 
surements. 



fraction measurement is important for the determination 
of the symmetric charge-breaking term D. In Ref. [3lj ]. 
a constrained fit to the measured branching fractions of 
J/tp — > BgBs is performed to extract the values of the 
parameters A, F, D' , F' and 5 using the Table [IV] am- 
plitude parametrizations. In the previous fit [3l|, D = 
was assumed, i.e., B{J/ip —> AS + c.c.) — 0. We per- 
form another fit that includes our new measurement and 
includes a nonzero value for D. The fit results are listed 
in Table [V] In comparison to the Ref. [3l| results, the 
value for the relative phase S has changed significantly, 
while the A, D 1 , F and F 1 values do not change signifi- 
cantly. The measurement of the isospin violating decay 
J/tp — > AS° + c.c. also provides useful information on the 
mechanisms for J/tp — > BgBio decays, where the large A- 
term is absent @, yj] . 

VI. SUMMARY 

In summary, with a sample of (225.2 ± 2.8) x 10 6 
J/tp events in the BESIII detector, the J/tp jAA 
decay has been studied. The branching fractions of 
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